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Abstract 


A charged particle array, named MATE-PA, which serves as an auxiliary detec- 
tor system to the Multi-purpose Active-target Time projection chamber for 
nuclear astrophysical and exotic beam Experiments (MATE) has been con- 
structed. The array is composed of twenty single-sided strip-silicon detectors, 
covering around 10% of the solid angle. It is dedicated for the detection of 
reaction-induced charged particles which penetrate the MATE active volume. 
The performance of MATE-PA has been experimentally studied using an alpha 
source, and a 36-MeV '4N beam injected into the chamber of MATE, filled with 
a mixture gas of 95% tHe and 5% COz under the pressure of 500 mbar, at 
the Radioactive Ion Beam Line in Lanzhou (RIBLL). The results demonstrate 
good separation of light charged particles with the forward double-layer silicon 
detectors of MATE-PA. The energy resolution of the Si detectors was deduced 
to be about 1% (ø) for an energy loss of about 10 MeV by the a particles. The 


inclusion of MATE-PA helps improve particle identification, and increases the 
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dynamic range for the kinetic energy of charged particles, in particular that of 
œQ particles up to about 15 MeV. 


Keywords: Silicon detector array, Active target, Time Projection Chamber 


1. Introduction 


Advances in the techniques for providing radioactive-isotope (RI) beams, 
and the constructions of high-intensity heavy-ion accelerator facilities 
Ml [5] worldwide have opened up vast opportunities to experimentally study and 
understand nuclei far removed from the beta-stability line [6] [7] [8] as well as 
nuclear reactions relevant to astrophysics. In an experiment that uses a reaction 
target, the selection of a target of appropriate material or phase (solid, liquid, 
or gas), as well as thickness is crucial to ensure luminosity, and to achieve the 
intended physical goals. Since the rates of RI beams are usually several orders of 
magnitude lower than those of stable beams, for nuclear spectroscopy in which 
charged particle(s) is(are) measured, it is essential to employ sufficiently thick 
reaction targets while attaining moderate energy resolution. Such requirements, 
however, are usually hard to be fulfilled with ordinary solid or liquid targets, 
since uncertainty in the reaction vertices and energy straggling will impair the 
resolution. 

The active-target time projection chamber (TPC) [9] [0]with a typical re- 
action gas thickness of 10 — 100 mg/cm? provides high luminosity necessary 
for nuclear spectroscopy with RI beams with little or moderate loss of resolu- 
tion. Since the reaction gas in an active-target TPC also serves as the detector 
medium, the reaction vertices of reaction-induced charged particles can be de- 
termined through tracking by the TPC. Because of these advantages, scores of 
active-target TPCs [I] [14] (15) 16, (17) (78) [79] have been constructed for 
experiments [20] [21] [22] [23] with RI beams. Recently, a prototype active-target 
TPC with 1000-channel readouts named Multi-purpose Active-target Time pro- 
jection chamber for nuclear astrophysical and exotic beam Experiments (MATE- 


1000) [24] was developed at the Institute of Modern Physics of the Chinese 
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Academy of Sciences (IMP-CAS), and applied to 7C+!*C fusion reaction mea- 
surement at stellar energies [24], and the energies above Coulomb barrier [25]. 
Building on the earlier success, a larger version of MATE with about 4000- 
channel readouts, denoted as MATE-4000 hereafter, has been constructed. The 
sensitive volume of MATE-4000 is 300(L)x300(W) x200(H) mm. To detect 
reaction-induced charged particles that penetrate the MATE active volume, a 
charged particle detector array, named MATE-PA (MATE Particle detector Ar- 
ray), consisting of twenty single-sided strip silicon detectors was designed and 
constructed. MATE-4000, as shown in Figfi{a), has been commissioned using a 
14N beam incident on *He, and employed in an experiment to measure inelastic 
alpha scattering of 11C to scrutinize the possible Z = 6 magic number [26] on 
the proton-rich side. 

In this paper, we report on the construction and performance test of MATE- 
PA. In Sec. [2] the design concept, geometry and readout electronics of the array 
are described. To assess the performance of the detector array, an experiment 
using '4N+“*He induced reactions has been performed. Monte Carlo simulations 
are also performed to estimate the geometrical acceptance of MATE-PA. Details 
of the experiment, as well as the experimental and simulation results are given 
in Sec. |3| The prospect of experiments with MATE-4000 in combination with 
MATE-PA are also presented. A conclusion and the future prospects are given 


in Sec. 


2. Charged Particle Detector Array for MATE 


In this section, the design concept and specifications of MATE-PA are pre- 
sented. The array was designed to achieve the largest possible angular coverage 
within MATE-4000. Detailed descriptions of MATE-4000 and its performance 


will be provided in another paper. 


2.1. Detector array 
The particle detector array consists of twenty single-sided strip silicon detec- 


tors, produced by Hamamatsu Photonics (Serial Nos. $10938-6734 and S10938- 
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Figure 1: (a) Side view of MATE-4000 with the constructed silicon detector array. The TPC 
field cage (sensitive area) is placed at the center of the chamber, with twenty silicon detectors 
mounted on the flanges around it. (b) Layout of the silicon detector array, seen from the beam 
direction. The side flanges have been tilted for clarity. See text for detailed descriptions of 


the assignments. 


6735), mounted in a single-layer 2 x 3 matrix on the two side walls, and a 
double-layer 2 x 2 matrix on the front wall of the target chamber. Each sili- 
con detector has eight strips with strip widths of 11.275 mm, strip gaps of 0.1 
mm, and a sensitive area of 91 x 91 mm/?; the size of one silicon is 100 x 100 
mm’. The side silicons have a nominal thickness of 600 um, while the front 
and rear forward silicons have thicknesses of 150 wm and 600 um, respectively. 
The upper (lower) rows of the right- and left-side silicon detectors are labeled 
RU; (RD;) and LU; (LD;), respectively, where the subscript 7 indicates 1, 2, 
and 3. The forward detectors consist of 2 layers of Si detectors, labeled as FF; 
and FS}, respectively, where the subscript j indicates 1, 2, 3, and 4 as shown in 
FigfI{b). For simplicity and to lower the risk of electric discharges, we mounted 
the silicon detectors close to the flanges. The distance from the RU; (RD;) and 
LU; (LD;) Si detectors to the field cage was set at 13.35 cm (14.00 cm), while 
that from FF; (FS;) was 17.80 cm (19.15 cm). The silicon detector array covers 
about 10.6% of the total (47) solid angle in the laboratory frame, assuming that 
the reaction point is at the center of the TPC. 

The purpose of choosing strip silicon detectors with moderately large strip 
widths, instead of large area single-channel or micro-strip silicon detectors, is 
to achieve sufficient energy resolution while keeping the readout channels at 


a manageable number. The current strip width is sufficient since the silicon 
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detector array will be used together with the TPC, which can provide track 
information, and the maximum rate of particles hitting the silicon detectors is 
expected to be of the order of a few hundreds per second. 

To protect the Si detectors from possible electric discharge, a layer of elec- 
trostatic shield (as shown in Figfif{a)) is placed between the MATE field cage 
and each side of the Si array. The shields were made by soldering Cu-Be wires 
of a diameter of 100 um at a constant distance of 1 cm on the copper-clad 
surface of a glass epoxy frame. Electrical grounding is achieved by connecting 
the shields to the vacuum chamber. Since no electric discharge was observed 
during the commissioning experiment, the distances will be optimized for future 
experiments to improve geometrical acceptance. For instance, the geometrical 
acceptance can be increased by about two times by reducing the distances be- 
tween the silicon detectors and the field cage from the current distance (around 


15 cm) to 5 cm. 


2.2. Electronics and data acquisition system 
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Figure 2: Typical electronic circuit diagram for the Si detectors. 


The output signals of the silicon detectors are fed to 16-channel integrated 
charge-sensitive preamplifiers, which were designed at China Institute of Atomic 
Energy. The preamplifier, named Smart PreAmplifier (SPA) [27], was designed 


as a compact and low-cost preamplifier with a fast response and a low power con- 
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sumption. Its smaller volume allows it to be placed closer to the detectors, thus 
reducing environmental noise, and enabling it to achieve similar performance to 
other commercial preamplifier modules. Due to the limited space in the MATE 
chamber, the SPA modules are mounted on the outside of the flanges. Every six- 
teen readout signals from a SPA module are fed to a Mesytec MSCF-16 module, 
which includes 16-channel shaping and timing filter amplifiers, as well as con- 
stant fraction discriminators, before being fed to an analog-to-digital converter 
(ADC; Mesytec MADC-32 module). The timing signals from the MSCF-16 
modules are sent to a time-to-digital converter (TDC; CAEN V1190A module). 
A typical electronic circuit diagram is given in Fig. B] The VME and GET (Gen- 
eral Electronics for TPC) [28] [29] based data acquisition system are employed 
for the data readout from the silicon and TPC detectors, respectively. Since the 
GET system acquires at a much slower rate (around 100 Hz for MATE), only 
the reaction events are selected to trigger the GET acquisition. In the current 
experiment, the OR signals from MSCF-16 modules and the GET system are fed 
to a gate generator (ORTEC GG8020) to provide the gate signals for VME and 
GET systems. Therefore, both DAQ devices start simultaneously and operate 


at the same rate. 


3. Experiment 


To evaluate the performance of the Si detector array, we have performed 
an experiment at the Radioactive Ion Beam Line in Lanzhou (RIBLL) [30] [31] 
at the Heavy Ion Research Facility in Lanzhou (HIRFL) [32] [33], operated by 
the IMP-CAS. A !4N primary beam with an energy of 117.6 MeV provided 
by the Sector Focusing Cyclotron (SFC) was transported to RIBLL. To reduce 
the energy of '4N so as to stop it in the TPC chamber, a 39.2-um- and a 
60.6-~m-thick Al plates were placed at the momentum-dispersive focal plane 
C1 and at the final achromatic focal plane T2, respectively. The resulted 14N 
beam with a total energy of about 36 MeV and an intensity of 10° pps was 


injected into the gas chamber of TPC, which was filled with a mixture gas of 
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95% tHe and 5% CO% at 500 mbar; the mixture gas acts as the target as well as 
the detection medium. Charged particles generated by reactions between 14N 
and the target nuclei were tracked and detected by the TPC, while part of the 
charged particles that penetrate the TPC were detected by the Si detectors. 
Data taking was performed using mainly the Si |J TPC trigger. To prevent 


the beam from triggering the TPC, the TPC sensitive area was divided into two 


regions: a low-gain region with a width of 1.7 cm around the beam axis for the 
detection of beam particles, and high-gain regions on both sides for detection 
of reacted particles and TPC trigger generation. 

For simplicity, the outputs of every four neighbouring channels of the Si 
detectors on the left and right sides were summed into one channel, which 
reduced the electronic readout channels without greatly affecting the energy 
resolution. Energy calibration for all Si detectors were conducted using a triple- 
nuclide alpha source, consisting of ?°9Pu (5.147 MeV), 244 Am (5.486 MeV) and 
244m (5.805 MeV), prior to the experiment. The typical energy resolution is 
about 0.8%(c) at around 5.5-MeV alpha particle energy. 


8.1. Particle identification 


3.1.1. Forward region 

The reaction-induced charged particles emitted in the forward direction were 
detected by the FF and FS silicon detectors. Figure Bla) shows a typical 
particle-identification plot obtained from the measured energy deposit in the 
first (AEfp) and second layer (AEps) of the forward silicon detectors. Three 
different loci corresponding to protons, deuterons and a particles are clearly 
seen. Here, we have eliminated overflow events and events due to accidental 
triggers taking into consideration the tracking information from the TPC. 

To better understand the data, simulations using the MATEROOT pro- 
gram [34] were performed assuming *He(!4N, p), *He(14N,d) and *He(14N, a) 
reactions, and taking into consideration the full geometry of the detector setup, 
the incident 14N beam energy, and the “He + CO2 mixture gas in the TPC 
chamber. MATEROOT is a FairRoot [35]-based data analysis and simulation 
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Figure 3: (a) AE-F plot obtained with the forward silicon detectors for charged particles 
induced by the *He+!4N reaction. The blue and green points represent experimental and 


(b) 


AE-AE plot after re-calibration, taking into consideration the penetration point of protons. 


simulated data, respectively. The axes are calibrated using a triple-nuclide a source. 


platform, built on CERN’s open source data analysis software ROOT [36] and 
the GEANT4 simulation framework [37], developed specifically for experiments 
with MATE. Results of the simulations (in green) are compared to the exper- 
imental data (in blue) in Fig. [3{a). Deviations between the experimental and 
simulation data are observed, which are likely to arise from insufficient energy 
calibration and inappropriate extrapolation along the energy-deposit axes. To 
eliminate the deviations, we consider the turnaround point of the proton locus, 
which correspond to protons that penetrate the second-layer silicons, in addi- 
tion to the triple-alpha-source data. The turnaround point was determined by 
searching for the crossing between two curves, obtained by fitting separately the 
loci for protons that penetrated and stopped in the second-layer silicons. The 
four data points, which include the turnaround point and the three data points 
from the triple-alpha source, were used to re-calibrate the energy-deposit axes. 
The resultant experimental loci after re-calibration are shown together with the 
simulation data in Fig. [3{b). The energy resolution of the Si detectors was 
deduced to be about 1% (c) for an energy loss of about 10 MeV by the a parti- 
cles. The results demonstrate sufficient particle identification for light-charged 
particles. To further discriminate protons from energetic deuterons that punch 


through both silicon layers, it is necessary to replace the second-layer silicon 
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by a thicker detector, e.g. CsI(Tl) or GAGG (Gadolinium Aluminium Gallium 
Garnet). 


8.1.2. Side region 

Most of the light-charged particles that penetrate the TPC and head towards 
the side silicons stopped in the single-layer silicon detectors under the current 
experimental condition. Figure [4] shows a typical scattering event recorded by 
the detector. The z axis in Fig. Ma) is defined as the beam axis, while the 
y axis is the axis perpendicular to the beam direction. z = 0 cm represents 
the upstream side of the field cage, or the entrance of the TPC sensitive region 
for the incident beam. The rectangles at the top and two sides of the figure 
represent strips of silicon detectors. The same event in the x(drift direction)-z 
plane is shown in Fig. ab). The recoil œ particle and the residual 14N particle 
are clearly observed simultaneously; the @ particle penetrates the TPC sensitive 
area, and hits the highlighted group of strips (in yellow color), while the 14N 
particle is stopped in the sensitive area. The charge and track information of 
the incident 14N beam particle for this reaction were not recorded because of 


the space charge effect at the beam region with low gain. 
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Figure 4: The trajectory of a typical scattering event on the (a) anode pad (y-z) plane, and 


(b) drift (x-z) plane. The z axis is the beam direction. 


For these events, particle identification (PID) can be performed using the 


charge (AErpc) and the track length (denoted as L) obtained with the TPC, as 
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well as the energy deposit in the side silicons (AEs;). The AErpc values were 
calibrated using the energy-loss curves simulated by the GEANT4 program. 
Typical AErpc-L and AErpc-AEsiprpc scatter plots are shown in Fig. 
For clearer identification of light and heavy particles, the AErpc axes of the 
PID plots are shown in logarithmic scale. The raw data plotted without the 
response of the front silicons are shown in Fig. la) and (b), in which two main 
components can be identified roughly as particles with atomic number Z > 
2, and light particles with Z < 2. It can be noticed from Fig. [B{a) that the 
expected light-charged particles with relatively high energy penetrate the TPC 
sensitive area from around 10 to 28 cm depending on the emitted angles, and 
some of them were detected by the side silicons, as shown in Fig. [5{b). 

The backgrounds can be attributed mainly to TPC electronic noise and 
possible electrons induced by the beam particles at the beam region around 
the center of the TPC sensitive area along the beam direction. The former 
can be eliminated by limiting the number of hits on the anode pad (e.g. 10 < 
Nnits < 200) of the TPC. The latter can be filtered out by applying geometrical 
constraints taking into account the kinematics of the two-body reaction, i.e. by 
requiring the positions of the first and last points of the trajectories to be within 
the high-gain region, and selecting appropriate track angles and directions for 
the two particles. A typical “good” event is shown in Fig.|4} After applying the 
above conditions, a clearer PID plot is obtained, as shown in Fig. [5{c) and (d). 
The backgrounds, especially the low A Erpc components, are removed, resulting 
in a clear separation between '4N and a particles. For the a particles that 
stopped in the TPC as shown in Fig. [{c), the energies down to a few hundred 
keV could be detected, which is consistent with the MATE-1000 result[24]. 
Thanks to MATE-PA, the detection energy range is extended to about 15 MeV 
for a particles as shown in Fig.[5{d). Protons and deuterons were barely detected 
due to insufficient gain of TPC in this experiment. A clump below the '4N band 


at around L = 20 cm in Fig.|5|(a) and (c) is due to a malfunctioning channel in 


the GET electronics, which resulted in quenched collected charges in a region 


close to the beam. 
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Figure 5: (a) Scattering plot of collected charge versus length of charged particles in the 
TPC, denoted as L. (b) AE-E plot for 14N and a particles measured with TPC (AE rpc) 
and silicon (AZsj).(c)(d) Similar plots after applying various conditions such as number of 


hits on the TPC anode pad, geometrical conditions, etc. For details, see text. 
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Figure 6: Correlation between the energy deposit in silicon and the angles of the particles 


obtained from Fig. Bo. 
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The results demonstrate feasible identification and separation of scattered 
light particles and residual heavy particles using MATE-PA incorporating the 
tracking information obtained with the TPC. It is worth mentioning that thanks 
to the tracking information, the angles of the particles can also be extracted as 
shown in Fig. a! which are useful for further analysis to refine the particle 
identification. 

For low energy charged particles that stopped in the TPC, one has to rely 
on the range information for particle identification. Another good alternative is 
to insert the TPC into a solenoid or dipole magnet, and identify particles based 


on their different curved paths in the magnetic field [14]. 


3.2. Geometrical acceptance 


For experiments that require measurements of absolute differential cross sec- 
tions, it is essential to determine the solid angle or geometrical acceptance 
of the detectors. To estimate the geometrical acceptance, simulations using 
MATEROOT were performed, taking into account the full geometrical setup, 
incident beam energy, medium gas, as well as the energy deposit of the beam 
and reaction-induced charged particles in the detectors. For thick-target ex- 
periments where the reaction energy changes drastically as the beam particles 
travel along the TPC, it is essential to iterate the geometrical acceptance sim- 
ulation taking into consideration the determined differential cross section. For 
simplicity, however, we assume an isotropic alpha scattering in the present work. 

For reactions within a very small range of excitation energy, geometrical 
acceptance can be represented as a function of polar angle in center-of-mass 
system. The simulated geometrical acceptance of the prototype MATE-PA for 
the 4N+*He elastic scattering is shown in Fig. [7{a). The geometrical acceptance 
e is defined as the ratio of silicon responses (total number of responses of 0 
or 1 per reaction event) to the number of reaction events generated within a 
subtended polar angle [0 — A0, 6 + A6], where 0 is the polar angle, and 2A6 
is its bin size. Only the events where TPC responded to recoil œ particles 


(L > 5 cm) were considered. Here, we have assumed a uniform distribution of 
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reaction vertices within the active volume of the TPC, and the reaction region 
is limited to the domain z = [0, 10] cm. For a more intuitive understanding, 
the geometrical acceptance and the solid angle (in unit steradian) are shown 
as functions of angle polar angle in the laboratory system in Fig. [7{b) by the 
left and right axes, respectively. The angular dependence of the geometrical 
acceptance arises mainly from the layout of the silicon detectors. 

To evaluate the dependence of the geometrical acceptance on the reaction 
vertex position, we divided the TPC sensitive area into ten parts along the z axis. 
For each part of the sensitive area, we obtained angle-integrated geometrical 
acceptance. The simulated result for the '4N+*He elastic scattering is shown 


in Fig. (7c), in which the shaded region is the TPC sensitive region. 
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Figure 7: (a) Geometrical acceptance € of MATE-PA for a particles induced by the 


4He(!4N,a)!4N reaction at different center-of-mass angles. The 14N beam energy at the 
entrance of the TPC chamber was about 36 MeV. (b) Geometrical acceptance (left axis) and 
solid angle in unit steradian (right axis) of the Si array at different laboratory angles. (c) Ge- 
ometrical acceptance for different reaction vertex positions along the beam axis; the shaded 


region is the TPC sensitive area. 


The average total geometrical acceptance is estimated to be about 6.0% for 
the *He(14N,a)!4N reaction. It is important to note that the energy of the 14N 
beam particles decreases quickly, and no reaction occurs at the end part of the 
sensitive area. This leads to a reduction in the overall geometrical acceptance. 
The actual geometrical acceptance, however, is expected to be larger. Based 
on the analysis of the experimental data, some of the reaction events occurring 
at the so-called “dead region” between the gas entrance window and the field 


cage were also recorded. Therefore, in the present simulation, we have also 
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expanded the reaction vertex position from the TPC sensitive area to the whole 
gas chamber. The results show that the reactions occurring near the front of 
the TPC sensitive area contribute a sizable geometrical acceptance. Hence, 
MATE and MATE-PA can be used to measure reactions that occur both within 
and outside the TPC sensitive area; moderate resolution might still be achieved 
with careful track reconstruction and proper track extrapolation. This aspect 


is very useful for thick-target experiments to measure excitation function using 


low-energy beam. 
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Figure 8: (a) Geometrical acceptance € of MATE-PA for a particles induced by the 


4He(12C,a)!?C reaction at different center-of-mass angles. The !?C beam energy at the 


entrance of the TPC was about 75 MeV/nucleon. (b) Geometrical acceptance (left axis) and 
solid angle in unit steradian (right axis) of the Si array at different laboratory angles. (c) Ge- 
ometrical acceptance for different reaction vertex positions along the beam axis; the shaded 


region is the TPC sensitive area. 


For comparison, a simulation was also performed for the elastic alpha scat- 
tering of a !?C beam with an energy of 75 MeV/nucleon incident on the TPC, 
filled with a mixture gas of 95% *He and 5% CO% at 500-mbar gas pressure. 
Results similar to those of Fig. [7] are shown in Fig. [S{a)-(d). The average total 
geometrical acceptance is estimated from Fig. [8{d) to be about 12.8%. 

The addition of MATE-PA increases the dynamic range for the kinetic en- 
ergy of charged particles, improves particle identification, and opens up broad 
opportunities for experiments measuring nucleon-transfer reactions, elastic and 


inelastic scatterings, charged-particle decay or evaporation. 
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4. Conclusion and future prospect 


A prototype particle detector array consisting of twenty single-sided strip sil- 
icon detectors for use with MATE-4000 to achieve large-angle and wide energy- 
range measurement has been constructed. The detector system was commis- 
sioned using a 36-MeV with a 10°-pps 14N beam incident on MATE-4000, filled 
with a 500-mbar mixture gas of 95% 4He and 5% COz at HIRFL-RIBLL. The 
test results demonstrate good separation of light charged particles with the for- 
ward double-layer silicon detectors of MATE-PA. The energy resolution of the 
Si detectors was deduced to be about 1% (øo) for an energy loss of about 10 
MeV by the a particles. Combining with MATE-PA, MATE-4000 will open up 
opportunities for experiments in nuclear astrophysics and RI beam physics in 
the near future. 

Recently, the readout cables of MATE-PA have been replaced by kapton 
ribbon cables to improve flexibility and work efficiency. The use of ribbon cables 
also helps improve environmental noise. In terms of the readout electronics, we 
use SPA and MSCF, incorporating a VME-based ADC for the data acquisition 
at present. We plan to migrate the General Electronics for TPCs (GET) system, 
which is used for the readout of the TPC signals, in the near future. An upgrade 
of MATE-PA with additional detectors such as CsI(Tl) or GAGG is also under 


consideration to accommodate a wide variety of experiments. 
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